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\^ ' Core-collapse supernova models depend on the details of the nuclear and weak interaction physics 

' inputs just as they depend on the details of the macroscopic physics (transport, hydrodynamics, 

■ etc.), numerical methods, and progenitors. We present preliminary results from our ongoing 

comparison studies of nuclear and weak interaction physics inputs to core collapse supernova 
, models using the spherically-symmetric, general relativistic, neutrino radiation hydrodynamics 

code Agile-Boltztran. We focus on comparisons of the effects of the nuclear EoS and the effects 
of improving the opacities, particularly neutrino-nucleon interactions. 
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1. Stellar core collapse and model inputs 

The evolution of massive stars (M > 10 M©) concludes when Si-burning builds a core of 
Fe-peak elements that can no longer be supported by electron degeneracy pressure and begins to 
collapse. During collapse, electron capture (EC) on nuclei and free nucleons reduces the electron 
fraction, Ye, in the Fe-core until the released neutrinos become trapped in the increasingly neutrino- 
opaque core. When the collapsing core reaches super-nuclear densities, ~ 3 x 10'^ g cm^^, 
nuclear repulsion halts further compression, and the core rebounds to form a shock, Mgh °^ Y}. The 
outgoing shock dissociates the heavy nuclei it passes, which, along with losses to neutrino emis- 
sion, saps the energy of the rebound shock, and it stalls above the newly-formed proto-neutron star 
(PNS). This standing accretion shock (SAS) is likely revived by heating from Vg and Vg absorption 
on free nucleons behind the shock, to disrupt the star as a supernova. Neutrino heating depends 
on the luminosity and thus on the opacities in the neutrino emitting region below. Multidimen- 
sional effects, including neutrino-driven convection and the SAS instability, expand the heating 
region and enhance the heating that drives the explosion. This paper is a preview of studies on 
the effects of changes in neutrino opacities [1] and the nuclear equation of state (EoS) [2] on the 
development of Fe-core collapse supemovae using Agile-Boltztran [3, 4, 1] — a spherically sym- 
metric, general relativistic, neutrino radiation hydrodynamics code. All of our models have 102 
adaptive mass shells; 20-energy-group, 8-angle-quadrature transport with updated opacities for all 
V species; LS 180 EoS [5]; and an up-to-date 15 progenitor [6] unless otherwise stated. 

2. Neutrino opacities 

The often used opacity set defined by Bruenn [7] includes emission, absorption, and scatter- 
ing on free nucleons and nuclei, scattering on electrons, and the pair source e+e^ o vv. Addi- 
tional interactions such as nucleon-nucleon bremsstrahlung, A^A^ o NNvv (where A'^ is any free nu- 
cleon) [8], and neutrino-pair flavor conversion, VaVa ^ VpVp [9], have been shown to be important 
sources of /XT-neutrinos. Replacing the independent particle model (IPM) for EC on nuclei with 
an emission table using detailed nuclear data [10, 11] (LMSH EC table) changes de-leptonization 
during coUapse and the p and Ye gradients that the shock traverses [12]. Improved interactions of 
neutrinos with free nucleons that include nucleon recoil, dense matter correlations, etc. [13, 14] 
and corrections for weak magnetism [15] reduce the v-nucleon opacities and provide for spectral 
change not available in the Bruenn [7] v-nucleon opacities. For these studies we have added to the 
opacity set from [7] the NN Bremsstrahlung, LMSH EC table, and enhanced v-nucleon opacities. 

We have implemented the Reddy et al. [14] v-nucleon opacities using an energy sub-grid 
without weak magnetism corrections and computed models with all permutations of the Reddy and 
Bruenn v-nucleon opacities. We find that the shock radius (Fig. [I]) is larger for models using the 
Reddy v-opacities (blue) than for the Bruenn v-opacities (black). The differential effect is larger 
for the change in scattering opacity than for the change in emission/absorption opacity by a factor 
of 3^ depending on the epoch [1]. The luminosities of Ve and Ve both increase with the smaller 
Reddy et al. [14] opacities, which allow neutrinos to stream more easily from the PNS while their 
r.m.s. energies remain largely unchanged as measured outside the shock. The difference between 
the increased luminosity from scattering and emission/absorption changes, respectively, parallels 
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Figure 1: Shock radii for 15 Mp, models [6] with different nucleon opacities plotted against time after core 
bounce. The blue line is the updated opacity set that includes the Reddy et al. [14] v-nucleon interactions. 
The black line uses the Bruenn [7] v-nucleon opacities. All other inputs and parameters are the same. The 
small jumps in shock radius result from the shock tracking algorithm shifting to an adjacent zone. 

the changes in shock radius, with the larger effect resulting from updated scattering. Both the 
increased luminosity and shock radii with updated nucleon interactions [13] were seen in the 1-D 
models of Rampp et al. [16]. 

3. Nuclear equation of state 

The nuclear EoS determines the thermodynamic properties of the fluid in SN simulations in 
two dynamically important regions — the high-density nuclear matter in the PNS and the nuclear 
statistical equilibrium (NSE) region at intermediate densities. Outside the Fe-core we include an 
optically thin non-NSE region with a Si/Fe mixture. We have computed simulations using two 
EoSs — the liquid-drop model of Lattimer and Swesty [5] (LS) EoS with a symmetry energy, E^y^ = 
29.3 IMeV and bulk incompressibility, = 180 and 375 MeV, and the relativistic mean field model 
of Shen et al. [17] (STOS) EoS with E^y^ = 36.9 MeV and = 281 MeV. A comparison of models 
using the LS EoS shows little variance with in initial shock position (~ 0.01 Mq, Fig. |2|, lower 
left) or core density (~ 5%, lower center) at bounce. At later times, the softer, more compressible, 
LS 180 EoS model with a more compact core had luminosities up to 10% larger (similar to the 
change from updated scattering) than the model with the stiffer LS 375 EoS ; and a small, detectable, 
change in the maximal extent of the shock (Fig. ^ upper). A similarly small difference due to in 
LS EoS models was seen by Swesty et al. [18]. 

The differences at bounce between the models using the LS and STOS EoSs originate pri- 
marily in the NSE region during collapse where the STOS EoS gives a lower abundance of heavy 
nuclei, resulting in less EC by nuclei [19] and therefore a larger initial shock position (0.49 Mq 
versus 0.44 Mq for the LS 180 EoS model. Fig. ^ lower). The differences in shock progression 
(Fig. ^ upper) are comparable to, but smaller than, those seen by [20, 19, 21], who used the LMSH 
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Figure 2: Upper panel: Shock radii for 15 Mq models [6] with different nuclear EoSs plotted against time 
after core bounce. Lower panel: Velocity (left), density (center), and Ye (right) profiles for the models at 
core bounce. The models use the Lattimer and Swesty [5] EoS with ffj = 180 MeV (blue, same as the blue 
line in Fig. 1) and fCj = 375 MeV (green), and the Shen et al. [17] EoS (red). 

EC as we did, but larger than the differences seen by [22] who did not and therefore Yg reduction 
during collapse was dominated by capture on free nucleons, the abundances of which are affected 
by £'sym- Each of the referred simulations used a different 15 Mq progenitor and a different opacity 
set which complicates the comparison. We also found that the differences for variations in EoS 
using the 12 Mq progenitor [6] are even smaller than they are for the 15 Mq progenitor [2]. 

4. Conclusions 

Improved v-nucleon interactions [13, 14] and EC on nuclei [11] show definite, and complex, 
effects on the evolution of the collapsed core and should be included in all realistic supernova sim- 
ulations. For example, we have shown here nontrivial and consistent enhancements in the shock 
radius by using the Reddy et al. [ 14] v-nucleon scattering opacities. Variations in the compressibil- 
ity of nuclear matter within the LS EoS seem to have little effect on shock progression through the 
times when multidimensional effects would start. Switching between the LS and STOS EoSs leads 
to larger differences than switching compressibilities within the LS EoS set, but these differences 
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remain smaller than those introduced by changing the neutrino opacities. Finally, the differences 
depend also on the progenitor used. Thus a systematic and comprehensive study must traverse the 
three axes of opacity, EoS, and stellar progenitor. Such a study is forthcoming [1,2]. 
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